Local Casimir energy for a wedge with a circular outer boundary 



Amir H Rezaeian 1 *and Aram A Saharian 2 t 

1 Department of Physics, University of Manchester 
Institute of Science and Technology (UMIST), 

PO Box 88, Manchester, M60 1QD, UK 

2 Department of Physics, Yerevan State University 
1 Alex Manoogian St, 375049 Yerevan, Armenia 

February 1, 2008 



Abstract 

The local Casimir energy is investigated for a wedge with and without a circular outer 
boundary due to the confinement of a massless scalar field with general curvature coupling 
parameter and satisfying the Dirichlet boundary conditions. Regularization procedure is 
carried out making use of a variant of the generalized Abel-Plana formula, previously estab- 
lished by one of the authors. The surface divergences in the vacuum expectation values of 
the energy density near the boundaries are considered. The corresponding results can be 
applied to the cosmic strings. 
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1 Introduction 



There are two main mechanisms which contribute to the energy of the vacuum: spontaneous 
symmetry breaking and the Casimir effect. The Casimir effect is a direct consequence of the 
quantum field theory and it is one of the most interesting macroscopic manifestations of the 
nontrivial properties of the physical vacuum. Although the Casimir effect has been extensively 
studied 0, [3|, |I| there are still difficulties in both its interpretation and renormalization 
[||, H . Moreover, the absence of a complete renormalization procedure, in practice, limits all 
calculations to the special case of highly-symmetric boundary configurations (parallel plates, 
sphere, cylinder) with a specific background metric. From this point of view the wedge with 
a cylindrical outer boundary is an interesting system, since the geometry is nontrivial and it 
includes two dynamical parameters, radius and angle, for phenomenological purposes. Due to 
the formal analogy that exists between a wedge and a straight cosmic string, the corresponding 
results can be applied to cosmic strings. Cosmic strings predicted in the framework of various 
gauge theories with spontaneously broken symmetries could have been created at cosmological 
phase transitions in the early universe [||, [7j. In the case of a static and straightline cosmic 
string the outside geometry is a locally flat conical spacetime with an angle 'deficit' 8tt/j> and /i 
is the linear mass density of the string. The superconducting strings form a particular subclass 
of cosmic strings [||. In the case of scalar field, it is reasonable to impose a Dirichlet boundary 
condition at r = a Q. The investigation of the influence of the cosmic strings on the behavior 
of quantized fields have been carried out in [1C]. In this paper we will study the Casimir energy 
density for a massless scalar field with a general curvature coupling parameter inside a wedge of 
opening angle cpo with and without a cylindrical outer boundary assuming the Dirichlet boundary 
conditions on the constraining surfaces. Some most relevant investigations to the present paper 



are contained in [11, [l2|, 13, 14] for a conformally coupled scalar and electromagnetic fields in a 



four dimensional spacetime. The total Casimir energy of a semi-circular infinite cylindrical shell 
with perfectly conducting walls is considered in [jll]] by using the zeta function technique. Our 
method here employs the mode summation and is based on a variant of the generalized Abel- 
Plana formula |l(| together with an adequate cutoff function. This allows to extract from the 
vacuum expectation value of the energy density the part due to a wedge without outer cylindrical 
shell and to present the cylindrical part in terms of the strongly convergent integrosum. We 
have organized the paper as follows. The next section is devoted to the consideration of the 
local Casimir energy for a massless scalar field with a general curvature coupling inside a wedge. 
This investigation generalizes the previously known result for a conformally coupled scalar and 
is essential for the main purpose of this paper in section ||, where the vacuum densities are 
considered for a wedge with the outer cylindrical shell. A formula for the shell contribution to 
the vacuum energy density is derived and the corresponding surface divergences are investigated. 
Finally, the results are re-mentioned and discussed in section ||. 



2 Vacuum energy density inside a wedge for a scalar field with 
general curvature coupling 

In this section we will consider a scalar field <p, with the curvature coupling parameter £, obeying 
Dirichlet boundary condition on the boundary of the wedge-shaped region formed by two plane 
boundaries intersecting at an arbitrary angle 4>o'- 

4>(t,r,4> = 0,2:1, ...,z N )= 4>{t,r,4) = 0o, Zi, ...,z N )=0, (1) 
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where we use cylindrical coordinates (r, (f>, Zi, . . . , zn) in D = N + 2 - dimensional space. The 
corresponding field equation is in form 



(V*Vi + £R) ip(x) = 0, 



(2) 



with R being the scalar curvature for the background spacetime. By using this equation, in the 
case of the flat background the corresponding metric energy-momentum tensor may be presented 
in the form 



£- \ ) .7//,V"V /( -iV,V, 



(3) 



The vacuum expectation values for these quantities can be derived by evaluating the mode sum 

(0\T ik (x)\0)=J2TikWa(x),<pUx)}, (4) 



where {(p a (x)} is a complete orthonormal set of positive frequency solutions to the field equation 
with quantum numbers a, satisfying the corresponding boundary conditions. In the region inside 
the wedge the corresponding eigenfunctions have the form 



1/2 



Ji{^r) sm(l<fi) exp Ukrii — iu;t) , 



a 



(2tt) n / 2 Wo. 
(n,7,k), < 7 < oo, — oo < kj < oo, n = l, 2, 



where rn = (z±, . . . , Zn), and Ji{z) is the Bessel function. Substituting the eigenfunctions into 
(|j) for the corresponding vacuum energy density one finds 



(Ou>|Too|0,j 



2(2vr) 



N, 



n=l 



£/ dr V I dNk l [ A i _) (7> k > r) - A\ +) ( 7 , k, r) cos(2/0) j (6) 



(+) 



with jO^,} being the amplitude for the corresponding vacuum inside a wedge, and we use the 
notations 
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A\ ±] ( 7 ,k,r)=^( 7 r)+(2£--) 7 Z 



1 



1 ± 



To evaluate expression (^) firstly we integrate over k by using the formula 



k s d N k 



1 



N+s-l 



7T 



(N-l)/2 



N + 8-1 



N + s 



^FTy ' r(JV/2) 

and next we integrate over 7 on the base of the formula 

.-a r (i^) r(/ + a/2) 



d7 7 a - i ^( 7 r) 



(7) 



(8) 



(9) 



20FT (1 - a/2) T(l + 1 - a/2) ' 
The integral involving the derivative of the Bessel function can be evaluated by using the relation 



1 



V 



sy 2 y«2 



d 2 1 d 
27 2 V cZr 2 r dr 



Jfiir). 



(10) 



As a result for the separate terms in (0) one finds 



/■«. /• —N/2 1 p f; , iV+i\ 

X ^y^k^( 7) k,r) = ^ r r(-iv/2-i) r (ii) 





[(Ar+l) 2 (/V + 2)(^-^)-Z 2 
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d-i I d N k^A\ +) (7,k,r) cos(2Z<£) 



7T 



JV/2-1 

— (N + 2)(Z-a c )T(-N/2 

2 



(12) 



(iV + 1) 



r(i 



cos(2/0), 



where £ c = (-D — 1)/(4D) is the curvature coupling parameter for the conformal case. To obtain 
the regularized value of the energy density we have to subtract from (Q) the corresponding 
quantity for the Minkowski vacuum without boundaries. The latter can be found in a similar 
way. The corresponding eigenfunctions are in form 



i 



7 \i/2 
2u) 



(2vr)(^+ 1 )/ 2 
Now on the base of formula (Q) 
(O a ./|T o(x)|Om) 



Jnirif) exp (in<t> + ikrii 



loot) 



n 



2, 



one finds 
1 



(27T) 



00 roo r 

iV+iE'/ *T /^k^H( 7 ,k,r), 

n=0 JO J ^ 



1,0,1,2.... 

(13) 
(14) 



where |0m) denotes the Minkowski vacuum, and the prime means that in the sum over n the 
term with n = is taken with weight 1/2. Taking into account that X^nLo '^n(l r ) = V 2 it can 
be easily seen that expression (|l4|) is presented in the standard form 

. r — I — l i \ I , /" d D \z oj 
{Om\T o(x)\Om) 



(2it) d 2 ' 



After integrating over k and 7 in (|i~4|) we receive 

rHv/2-i) ~ r(»+^; 



(Om|T o(x)|Om) 



2JV+3 7r JV/2+2 r iV+3 



E- 

^ r 



iV-l 



[(7V + l) 2 (iV + 2)(C-^) 



n 



(15) 



(16) 



n=0 "V" 2 1 

Let us consider the case N = 1. The sums over n are evaluated by using the Riemann zeta 
function ((s) = n~ s , C(-l) = -1/12, C(-3) = 1/120, and the formula 



n p cos(/9ra) 



n=l 



1 irp d p , n 
-sm— — — cot(/3/2 , 

2 2 d/3P VM/ ; ' 



(17) 



valid for an odd p. For the regularized vacuum energy density this leads 



00/rcg 



(O w \Too\O w ) 



+ 



(0m\T O0 \0 m ) 

_2 



1 



csc (tt 



■0 



7T 



2cpl 



1440vr 2 r 4 
3csc 2 (7T</>/0 o ) - 2) + 1 



(tt/^o 



12vr 2 r 4 



+ 



(18) 



where £ c = 1/6 is the curvature coupling parameter for the conformal case in D = 3. In formula 
(18) the first two terms on the right come from the summand with ^ ( 7 ,k, r) in (^), and 
the last term comes from the summand with A^ ( 7 ,k, r). For <pQ = it it can be easily seen 



that expression ( ^L8|) coincides with the well-known result for the Casimir energy density in the 
case of a single plate geometry. For the conformally coupled case from (^) we obtain the 
standard result previously derived in Refs. |12|,[fT8|. It is of interest to note that in this case 
the renormalized energy density is angle independent and is finite on the sides of the wedge. 
It diverges on the edge as r -4 , when r — > (as r~( D+1 ^ for a D - dimensional space). For 
a non-conformally coupled scalar the energy density is angle dependent and contains surface 
divergences on the wedge sides. Note that for a minimally coupled scalar (£ = 0) the regularized 
energy density (|i~8|) is negative everywhere, (2oo) re g|g=o < 0. 
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3 Vacuum energy density inside a wedge with circular boundary 

Having investigated the vacuum energy density inside a wedge we now turn to the case with 
additional circular boundary with radius a. The boundary conditions are Dirichlet ones: 



<p(t,r,<p = 0,z 1 , ...,zn) = <p(t,r, 
ip(t,r = a,4>, zi, . . . ,zn) = 0. 



f) ,zx, ...,z N ) = 0, 



(19) 



The eigenfunctions satisfying these boundary conditions are of the form 

(faix) = fiaJidr) sin(Z</>) exp («krii — iujt) , 



(20) 



where the notations are the same as in (|5]) . The normalization coefficient (5 a is determined from 
the standard Klein-Gordon scalar product and is equal to 



0i 



(2vr) Jv w(/ ) oa 2 J/ 2 (7a)' 



(21) 



The eigenvalues for the quantum number 7 are quantized by the boundary condition (19) on 
the cylinder surface r = a. From this condition it follows that the possible values of 7 are equal 
to 



7 = A n ,-/o, J = 1, 2, • • • , 



(22) 



where X n j are positive zeros of the Bessel function, Ji(X n ,j) = 0, arranged in ascending or- 
der, X n j < A n j + i. Substituting the eigenfunctions (|20|) into mode sum formula @ for the 
corresponding energy density one finds 



n=l j=l 



(0|r 00 |0) = ~ ^ I A V (W«. k ' r ) " ( W«> k, r) cos(2^) 



(+) 



(23) 



with a set of quantum numbers a = (n,j, k). The vacuum expectation value ( p3[ ) is divergent. 
To make it finite we introduce the cutoff function tp^d) satisfying the condition tp^ — > 1, jj, — * 0. 
To extract the divergent part we will apply to the sum over j the summation formula |l6| 



/(Anj 



^n,jJ'i 2 (K,j)J A 2 j + C 2 



7T, 



--dz + —Res 



z=0 



MYi(z) 
Vz 2 + c 2 Ji{x 



+ 



2tt 



dz 



Ki(z) e- lni f(ze ni / 2 ) + e lvri f{ze~ ni / 2 ) 
Il(z) 

Ki(z) e~ l ' Ki f{ze Ki l 2 ) - e M f{ze" Ki / 2 ) 



\/ z 2 — c 2 



, (24) 



where Yi{z) is the Neumann function, and Ii(z), K[(z) are Bessel modified functions. This 
formula is valid for functions f(z) satisfying the conditions 



\f(z)\ <e(x)e c ^\, z = x + iy, a < 2, 



/(*) 



o\z 



0, 



(25) 
(26) 



where e(x) — ► for x — ► 00. 
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To evaluate the sum over j in (23) now as a function f(z) we choose 

f{z) = z \a^ ] (z/a, k, r) - 4 +) 0/a, k, r) cos(2/0)j %{z/a), 



(27) 



assuming a class of cutoff functions for which f(z) satisfies conditions ( |25D and fl26|) uniformly 
with respect to the cutoff parameter [i. It can be seen that the contribution of the first integral 
on the right of formula (|24| ) to the energy density (^) corresponds to the energy density for 
the case of a wedge without outer cylindrical boundary. This quantity was investigated in the 
previous section. The second and third integrals on the right-hand side of fl24l ) are finite in the 
limit fj, — ► 0. Using the standard formulae for the Bessel function we see that the subintegrand of 
the second integral is proportional to ip^iz/a) — ip^—iz/a). Consequently, after removing the 
cutoff the contribution of this integral will be zero. Hence, omitting this integral and removing 
the cutoff for the vacuum energy density we obtain 



(0|T 00 |0) 



(OjToolO 

V 



1 



2 N ~ 1 7T N + 1 C 



n=l 



d N k 



zdz 



Ki(z) 



ka Vz 2 - k 2 a 2 h{z) 



k 2 I 2 (zx) 
2C 



1\ l 2 lf(zx) 



a 2 x 2 



cos (l(f>) > , x 



r 
a 



2g + l/2 
2^-1/2 



sin 2 (14>)- 



(28) 



The integration over k can be done using the formula [19] 



k s g(z)dz 



7T 



N/2 



k v^2TTfc2 r(JV/2) 



B 



N + s 1 
2 '2 



dzz 



N+s-l 



9(z), 



(29) 



where B(x,y) is the Euler beta function. As a result for the vacuum energy density (|28| ) one 
obtains 

(0|T 00 |0) = (0 W \T 00 \0 W ) + (Too)c, (30) 

where (0^,1 Too |0 W ) is the vacuum energy density for a wedge without outer cylindrical boundary, 
and the term 



4^ - 1 



00/c 



2 D-l 7T (D-l)/2^ a D + l T{{D _ 1)/2) Ey 



dz z 



D 



Ki(z) 



B\ ) (zx) - B^> (zx) cos(2Z0) 

(31) 

is the contribution due to the presence of the cylindrical shell at r = a. Here I is defined in 
accordance with (IBI) and we have introduced the notations 



(+) 



B\ ±] {y) = I[ 2 (y)+lf(y) 



2 f_ 

1 + (D-im-i) T ^ 



(32) 



In accordance with the problem symmetry, expression ( |3l| ) is invariant under the replacement 
(j) — > 4>q — (p. Aiming to compare with the result for the energy density of a cylindrical shell with 
the radius a let us write down the corresponding formula, which is obtained from the general 
result of [19] and has the form 



(^bo)reg " — 



2 ^ 7r (D+i)/2 a D+i r (( jD _ !)/ 2 ) 



E 



iIzz D^$ B ( r ) {zx) 



In{z) 



(33) 
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with the same notations as in (|3l|). 

For < r < a the cylindrical part ( |3l"| ) is finite for all values < <fi < <fro, including the 
wedge sides. The divergences on these sides are included in the first term on the right-hand side 
of (pC), corresponding to the case without circular boundary. At the edge r = the boundary 
(3l|) vanishes for < <pQ < tt, is equal to 



part 



4£-l 



00/0|r=0 



2 D 7r (D+l)/2 r (( jD _ l)/ 2 ) a B+l J Q 



dzz 



D 



Ki(z) 
Hz) 



TT 



(34) 



for cf)Q = tt, and diverges as r~ 2 ( 1_7r /^°) for 0o > tt. 

The boundary part {Tqq) c diverges on the cylindrical surface r = a as well. To investigate 
the corresponding behavior near this surface let us consider the integrosum in (|3l"l), introducing 
a new integration variable z — * lz: 



I= Y^i D+i 



71=1 



dz z 



D 



Ki{lz) 
k{lz) 



B\ ) (Izx) - B\^> (Izx) cos(2/0) 



(+) 



(35) 



By taking into account that near the surface r = a main contribution comes from the large 
values of I we can replace the Bessel modified functions by their uniform asymptotic expansions 
for large values of the order (see, for instance p7j| ). In the leading order over 1/7 one finds 



1 



where 



1 

2^2 



n=l 



D 



dz z D ~ 2 (l + z 2 x 



2 rr 2 ) l/2 e -2l[ V (z)- V (zx)] 



b { ~\zx) - b {+) (zx) cos(2Z0)l , (36) 



b {±) (y) 



1 + 



2y 2 



(Z>-1)(4C-1) 



rj(z) = yl + z 2 + In • 



l + VTTz 2 

Near the cylindrical boundary in the leading order over I — x from (^) we obtain 

D oo 



(37) 
(38) 



1 / TT 



E 

n=l 



J? 



D 



dz Z 



D-2 



2\l/2 -2n7r(l-x) v / T+l2/ fa 



(l + z'Y^e 



o 



- 6 (+) (z) cos(2Z0) 



On the wedge sides cos(2Z(/>) = 1 and this yields 

rco dzz D ~ 2 
:"'./, OT^ 6 

Summing over n by using 

<i D 1 



(\ D co 



2n7r(l-x')v / l+2 2 /^o 



0, 



(39) 
(40) 



E 

n=l 



„f —an 

n e 



(-1) 



D 



T{D + l) 



da D 



a 



D+l 



a«l 



(41) 



we have 



(f> r(D + 1) 



D-2 



M D-1)T 2 (^i) 



2^+i 7 r(l-x) c + 1 y (l + z 2^/2+i" I6vr(l-x; 



D+l 



0, = 0o- (42) 
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Substituting this into formula ([?!]) to the leading order over ( 

(4£-i)rm 



L one finds 



00/c 



2 D+2 7r (Z>+l)/2( a _ r ) 



D+l ' 



0, 



(43) 



For D = 3 and a conformally coupled scalar field this term coincides with the first summand on 
the right of formula ( |l8|) with the distance (a — r) from the edge r = a, cj) = (or = 0q) and 
the opening angle ir/2. 

For the angles < < <po by using the formula 



oo 

£ 

71=1 



cos n/3 



(-1) D d D 



sinh a 



2 da D \ cosh a — cos (3 



(44) 



introducing a new integration variable y = 2n(l — x)VT + z 2 /4>q and expanding over (1 — x) one 
finds that the leading contribution of the term with &( + ) (zx) cos 210 into (|36| ) is equal to 



2^+1(1 



1 + 



(D-i)(4e-i) 



(45) 



In this case the leading contribution from the term with i>( \zx) dominates and to the leading 
order by the way similar to that for ( |43"|) one has 



00/c 



g(H)r(f) 

2 D 7r (D+l)/2( a _ r )D+l' 



< < 



(46) 



This leading divergence coincides with the corresponding one for a cylindrical surface of the 
radius a (see, for instance, ]i~9|]). 

The surface divergences in the renormalised vacuum expectation values of the local physical 
observables result from the idealization of the boundaries as perfectly smooth surfaces which are 
perfect reflectors at all frequencies, and are well known in quantum field theory with boundaries. 
They are investigated in detail for various types of fields and general shape of smooth boundary 
18|. Near the smooth boundary the leading divergence varies as (D + l)th power of the 
distance from the boundary. For conformally invariant fields the coefficient of this leading term 
is zero. For non-smooth boundaries such as here, the latter is not the case and this leads to the 



extra divergences in the global quantities such as total Casimir energy (see, for instance, [15| for 
the case of a semi-circular infinite cylinder and discussion in p0| ). It seems plausible that such 
effects as surface roughness, or the microstructure of the boundary on small scales (the atomic 
nature of matter for the case of the electromagnetic field pi] ) can introduce a physical cutoff 
needed to produce finite values of surface quantities. An alternative mechanism for introducing 
a cutoff which removes singular behavior on boundaries is to allow the position of the boundary 
to undergo quantum fluctuations [^2| . Such fluctuations smear out the contribution of the high 
frequency modes without the need to introduce an explicit high frequency cutoff. 

The dependence of the cylindrical boundary part of the vacuum energy density multiplied 
by a D+1 , a D+1 

(Too) a ° n the angle 4> and radial coordinate x = r/a is shown in figure [y for a 
minimally coupled scalar (£ = 0) and in figure || for a conformally coupled scalar (£ = £ c ) in 
D = 3. The cases (po = ir/2 (left graphics) and <f>o = vr/4 (right graphics) are presented, when 
the cylindrical part vanishes at the edge r = 0. 
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Figure 1: The expectation values for the cylindrical part of the vacuum energy density, 
a D+1 (Too) c for D = 3 minimally coupled scalar field in the cases 0q = it/2 (left) and cfto = 7r/4 
(right), versus (ft and x = r/a. 




Figure 2: The same as in fig. |] for a conformally coupled scalar. 
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4 Final remarks 



In the present paper we calculate the vacuum energy density for a massless scalar field with a 
general curvature coupling and satisfying the Dirichlet boundary conditions in the wedge-shaped 
region formed by two plane boundaries intersecting at an arbitrary angle 0o with and without 
a cylindrical outer boundary. All calculations are made at zero temperature and it is assumed 
that the boundary conditions are frequency independent. The latter means no dispersive effect 
is taken into account. The energy density for the wedge, given by formula (|i~8|), has been 
obtained by standard application of the mode summation method. This formula generalizes the 
result previously known in literature for a conformally invariant scalar field. In the case of a 
non-conformally coupled scalar the vacuum energy density is angle-dependent and diverges on 
the wedge sides. For a minimally coupled scalar it is negative everywhere inside a wedge. In 
the case of presence of an outer cylindrical boundary, discussed in section || the expectation 
value of the corresponding vacuum energy density is presented in the form of series over zeros 
of the Bessel function. The summation formula for this type of series based on the generalized 
Abel-Plana formula allows to extract explicitly from the expectation value the part due to the 
wedge without a cylindrical shell. The additional cylindrical contribution to the vacuum energy 
density is presented in the form of the strongly convergent integrosum. The latter is finite on 
the wedge sides but contains surface divergences on the cylindrical boundary. At the edge r = 
it vanishes for < 4>q < n, is finite for <po = ir, and diverges as r - 2 ( 1 - 7r /? i o) f or ^ > 7r _ The 
generalization of the results obtained here for the Neumann, or more general Robin boundary 
conditions is straightforward. For instance, for the Neumann case in the expressions (||) and ( p0|) 
of the eigenfunctions the function cos(7(/>) stands instead of sin(/(/>) and the quantum number n 
takes the values 0, 1, 2, . . .. In the case with a circular outer boundary now the eigenvalues for 
7a are zeros for the derivative of the Bessel function. The formula to sum the series over these 
zeros can be taken from fl6[| . 
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